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Abstract 

We study the interplay between Dirac eigenmodes and center vortices 
in SU(2) lattice gauge theory. In particular we focus on vortex-removed 
configurations and compare them to an ensemble of configurations with 
random changes of the link variables. We show that removing the vortices 
destroys all zero modes and the near zero modes are no longer coupled 
to topological structures. The Dirac spectrum for vortex-removed con- 
figurations in many respects resembles a free spectrum thus leading to a 
vanishing chiral condensate. Configurations with random changes leave the 
topological features of the Dirac eigensystem intact. We finally show that 
smooth center vortex configurations give rise to zero modes and topological 
near zero modes. 
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1 Introduction and outline of results 



The QCD vacuum has a highly non-trivial structure. Its excitations give rise 
to the key features of QCD, namely confinement and chiral symmetry breaking. 
Over the years, for both features several mechanisms have been proposed. For 
confinement the two most influential pictures are based on monopoles, respec- 
tively center vortices, while for chiral symmetry breaking instanton-type excita- 
tions have played a major role. It is widely expected that confinement and chiral 
symmetry breaking are linked through some unifying mechanism. This belief is 
supported by numerical results indicating that at the QCD phase transition de- 
confinement and restoration of chiral symmetry take place at the same critical 
temperature. However, our understanding of such a unifying mechanism is still 
in its infancy. 

Lattice QCD is an important tool for analyzing excitations of the QCD vac- 
uum and is even necessary for formulating some of the mechanisms such as the 
picture based on center vortices. Lattice simulations provide one with thermal- 
ized configurations of the gauge field as they appear in the path integral. These 
configurations can then be analyzed and one can try to identify the excitations 
of the gluon field relevant for different physical features. 

Such lattice investigations have provided strong evidence for the center vortex 
picture of confinement. Vortices obtained after center projection in the maximum 
center gauge PP, are physical (in the sense of the renormalization group) j21, and 
seem to constitute the relevant infrared degrees of freedom of Yang-Mills the- 
ory. Indeed, when center vortices are removed from the Yang-Mills ensemble jHI , 
the confinement properties are lost |3l E] . The loss of confinement is indicated 
by the vanishing of the string tension and the change of the infrared behav- 
ior of the Green's functions in Landau gauge E] and in Coulomb gauge [7j. 
The emergence of a string tension can be easily understood in the center vor- 
tex picture. In addition, center vortices provide also an appealing picture of the 
finite temperature deconfinement phase transition as a depercolation transition 
in a 3-dimensional slice of the lattice universe, taken at a fixed spatial coordi- 
nate |H1 El- Moreover, center vortices reproduce not only the correct order of 
the phase transition but also its universality class [TU]. On the other hand in 
the spatial 3-volume, there is no depercolation transition and the center vortex 
ensemble correctly describes the increase of the spatial string tension above the 
critical temperature fTTl. 

All these properties of center vortices detected on the lattice after center 
projection, as well as the correct order of the deconfinement phase transition 
are well reproduced in a random center vortex model for both SU(2) and 
SU(3) ^3]. Furthermore, the vortex percolation in the confinement phase is 
consistent with the vanishing of the free energy of center vortices in this phase 
as observed on the lattice Jl] and in the continuum |15j . 

When center vortices are removed from the Yang-Mills ensemble, not only the 
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string tension is gone but also the quark condensate vanishes and the topological 
charge is lost OE]. These findings suggest that center vortices could be respon- 
sible for all infrared key features, confinement, chiral symmetry breaking and the 
chiral anomaly. 

In D = 4 center vortices are closed fiux surfaces, which act as wave guides 
for low-energetic quarks ^H] and their percolation very likely triggers the con- 
densation of quarks. Their topological charge can be understood in terms of 
their intersection number flEl [W\ or, when center vortices are considered as time- 
dependent closed loops in 3-dimensional space, by their writhing number 

For the study of long range structures and topological excitations of the QCD 
vacuum, quantum fiuctuations pose a serious challenge in the analysis of lattice 
configurations. These hard UV modes dominate the action and the long range 
structures are hidden under short distance noise. In recent years it has been 
understood that the low-lying eigenmodes of the lattice Dirac operator provide 
a natural filter sensitive to long range structures. These low-lying modes can be 
studied and information on the underlying IR structures of the gauge field can 
be extracted. 

In this article we study the connection between center vortices and properties 
of Dirac eigenmodes for quenched SU(2) configurations. In particular we apply 
the technique of Ref. P] to remove the center vortices. 

A central question of our article is how the removal of center vortices affects 
the long range topological structures of the Yang-Mills vacuum. As outlined, this 
question can be addressed through an analysis of the low-lying Dirac eigenmodes. 
We will show in this article that topological modes are destroyed when removing 
the center vortices. In particular the zero modes are gone completely. When 
analyzing the low-lying modes with non-vanishing eigenvalues, we find that their 
local chirality is gone and they do not resemble small perturbations of instanton- 
type zero modes. 

For the spectrum of the Dirac operator we show that in many respects the 
spectrum of vortex-removed configurations resembles the spectrum of free fermions 
This implies that the vortex-removed configurations cannot build up a non- 
vanishing chiral condensate via the Banks-Casher relation. 

A legitimate criticism of removing the vortices is the fact that this procedure 
is a quite drastic modification of the gauge field. In this article we implement 
a crucial test by applying random changes to the original gauge configurations. 
In particular we multiply randomly chosen links of the lattice with the non- 
trivial center element. We demonstrate that to a large extent the topological 
information is stable under such random changes. The number of zero modes 
is essentially invariant and also the chiral properties of zero-modes and near 
zero-modes. Our test shows that the center vortices are correlated in a highly 
non-trivial way and their removal specifically destroys topological features of the 
gauge field configuration. 

Finally we will demonstrate for smooth vortex configurations that they indeed 
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give rise to zero modes of the Dirac operator. We analyze the procedures of 
center projection and vortex removal for these configurations. The findings for 
the smooth configurations support our interpretation of the results for center 
projection and vortex removal applied to thermalized configurations. 

Our article is organized as follows: In the next section we discuss technical 
aspects of our calculation, in particular the preparation of the original, the vortex- 
removed and the random-changed configurations, as well as the Dirac operator 
and the computation of the eigensystem. In Section 3 we discuss eigenvalue 
spectra and their interpretation. In Section 4 we present the observables based 
on the Dirac eigenvectors. In Section 5 we study Dirac spectrum and eigenmodes 
for the hand-constructed smooth vortex configurations. Our paper closes with a 
summary and the discussion of the results. 

2 Technicalities 

2.1 Preparation of the original ensembles 

The SU(2) lattice configurations were generated with the standard heat bath 
algorithm using the Wilson action. After a careful thermalization, each "mea- 
surement" was taken after 20 dummy heat bath sweeps. We analyze a total of 
100 quenched SU(2) configurations generated on a 12^ lattice at ^ = 2.5. 

2.2 Center vortices and their removal 

The method of identification IJ and removal of center vortices is based on the 
so-called maximal center gauge (MCG) [H [T7] . If 



denotes the gauge transformed link, MCG fixing is implemented by maximizing 



with respect to Q{x) thereby bringing each link as close as possible to a center 
element jEj, or a given gauge field configuration as close as possible to a collec- 
tion of center vortices PH]- This gauge condition was implemented by using an 
iteration-overrelaxation (10) procedure (details are presented in jl7]). The 10 
procedure was stopped when the difference between the variational action Sa^ of 
two subsequent gauge fixing sweeps was smaller than 10~^. 

After MCG fixing the center vortices are identified by replacing each link 
U^{x) by its closest center element ^^(x) G Z{2), thereby each gauge configura- 
tion U^{x) is converted into an "ideal" center vortex configuration, consisting of 
closed surfaces of plaquettes being equal to a non-trivial center element 1. 
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The center vortex-removed theory [Hj is defined by replacing the original gauge 
fixed links U^{x) by 

uHix) - z;{x)u;!{x) , (3) 

where the Z^{x) are the center projected counter parts of the U^{x). Let us 
emphasize that the vortex removal procedure [Sj de facto removes the center 
projected image Z^j_{x) from the original gauge fixed configuration U^{x). 

2.3 Configurations witii random changes 

One has to keep in mind that removing the vortices is a drastic alteration of 
the gauge configuration: A large portion of the links is fiipped, i.e. the links 
are multiplied by a factor of —1. This change leads to a fiip of about 3% of 
the plaquettes. The average plaquette drops by about 5 %, from 0.651 on the 
original configuration down to 0.621 after removing the vortices. It must be 
stressed, that the non-trivial center elements are highly correlated beyond the 
fact that the vortices form closed surfaces on the dual lattice. An important 
question is what the effects of a similar, but unstructured alteration of the lattice 
configuration will be. 

To address this question we analyze the effect of uncorrelated random changes. 
The procedure for these random changes is somewhat similar to the prescription 
for removing the vortices. We first randomly choose a certain number NfHp of 
links {x, /x) where the random changes are implemented. For the selected links 
we change the corresponding link variables U^{x) according to 

U,{x) ^ -U,{x) , (4) 

i.e., we fiip a randomly selected set of NfHp link variables. The number NfUp of 
randomly chosen links subject to the transformation Q is a free parameter. 

The number of fiipped plaquettes increases quickly with the number Njnp of 
randomly fiipped links. We find that randomly fiipping NfUp = 1000 links on 
our 12^ lattice (this is 1.2% of the total number of links) leads to a fiip of about 
4.6% of the plaquettes. This decreases the average plaquette by about 10%, from 
0.651 on the original configuration to 0.590 after the random changes. 

We remark that a simple combinatorial calculation gives the fraction p of 
fiipped plaquettes when fiipping a fraction q of links as p = 4g — 12g^ + 0{q^). 
Setting q = 0.012 (= 1.2 %), the value for our random changed ensemble, one 
finds p = 0.046 (= 4.6 %), as we observe. Since the links are chosen randomly, 
in leading order the fiipped plaquettes are unbiased. Thus we can estimate the 
average plaquette after the random changes to be (1 — 2p)0.651 = 0.591, almost 
the number we observe (0.590). If on the other hand one inserts p = 0.03 (= 3 
%), the value for vortex removal, one expects an average plaquette of 0.612. The 
discrepancy of this number to the observed value of 0.621 for vortex removed 
configurations, is due to a non-trivial correlation of the vortices. 
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2.4 The eigensystem of the Dirac operator 

For all four ensembles, original, center projected, vortex-removed and random- 
changed, we computed eigenvalues and eigenvectors of the chirally improved (CI) 
lattice Dirac operator [201 E]- The CI operator is a systematic approximation 
of a solution of the Ginsparg- Wilson equation which governs chiral symmetry 
on the lattice. The approximation is ultra-local and so allows for good chiral 
properties at a relatively low cost. The ultra-locality makes the CI operator also 
particularly suitable for the analysis of topological objects. In j22I it was e.g. 
shown that the CI operator provides a better lattice image of the continuum zero 
mode of instantons small in lattice units, when compared to the exactly chiral, 
but non-ultra-local overlap operator. 

The actual calculation of the eigensystem was done with the implicitly restarted 
Arnoldi method jSHI- For each configuration we calculated the 50 smallest eigen- 
values and the corresponding eigenvectors. We vary the boundary conditions of 
the Dirac operator which is a powerful tool to probe the system |21I-|2Z|- While 
we keep periodic boundary conditions for the spatial directions, we allow for an 
arbitrary phase in the temporal direction. Thus the eigenvectors v obey 

v{x + Li) = v{x),i=l,2,3 ; v{x + LZ) = e'^'''^ v{x) , (5) 

where L is the size of our lattice and i, 4 are the unit- vectors in the spatial, respec- 
tively the temporal directions. The phase factor ip can assume values between 

and 1, and the special cases of (/9 = and ip = 1/2 correspond to periodic, 
respectively anti-periodic temporal boundary conditions. 

For a Dirac operator which obeys the Ginsparg- Wilson equation, the eigen- 
values are restricted to the so-called Ginsparg- Wilson circle, a circle with radius 

1 and center (1, 0) in the complex plane. For the CI operator, which is an approx- 
imate Ginsparg- Wilson operator, the eigenvalues fluctuate around the Ginsparg- 
Wilson circle. In particular also the zero modes do not have eigenvalues that 
are exactly zero. In general they have a small, but non- vanishing real part. We 
stress, however, that there is no mixing of zero modes and near zero modes, since 
eigenvectors that correspond to eigenvalues with non-vanishing imaginary parts 
have an exactly vanishing matrix element with 75 and can be identified unam- 
biguously. Thus eigenmodes with real eigenvalues are the would-be zero modes 
of the continuum. 

3 Spectra of the Dirac operator 

3.1 Spectra for center-projected configurations 

As a first approach to analyzing the relation between center vortices and Dirac 
eigenmodes one can calculate the spectrum directly for the center projected con- 
figurations. The result of such a calculation is shown in Fig. ^ where in the 
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Figure 1: Dirac eigenvalues X in the complex plane. We compare the spectrum 
for the original ensemble {Ihs. plot) to the spectrum for center-projected conhgu- 
rations {rhs. plot). Each plot contains the 50 smallest eigenvalues for 10 different 
configurations with every configuration represented by a different symbol. Anti- 
periodic temporal boundary conditions were used for the Dirac operator. 



right-hand side (rhs.) plot we superimpose spectra for 10 center projected con- 
figurations and compare them to the corresponding spectra of the original con- 
figurations (Ihs. plot). All spectra were calculated using anti-periodic boundary 
conditions. 

In the Ihs. plot for the original configurations one clearly sees that the eigen- 
values cluster near the Ginsparg- Wilson circle (represented by the dashed curve) 
and only slightly scatter away from it (except for some real modes). The density 
of eigenvalues extends all the way to the origin indicating, that chiral symme- 
try is broken according to the Banks-Casher formula j2H]- In addition we find 
several real eigenvalues, i.e. would-be zero modes corresponding to non-vanishing 
topological charge via the index theorem. 

The spectra for the center projected configurations on the rhs. plot show a 
completely different picture. The spectra have developed a large gap, and all 
real eigenvalues are gone. The eigenvalues are concentrated in two clusters (sym- 
metric with respect to refiection on the real axis). These clusters are shifted 
to large values of the imaginary parts, beyond the range for the eigenvalues for 
the original configurations. This hard to interpret outcome is not really a sur- 
prise: The Dirac operator, which contains the gradient operator, can detect only 
smooth topologically non-trivial structures. The center projected configurations 
are, however, maximally discontinuous. The link variables can only jump from 
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Figure 2: Dirac eigenvalues A in the complex plane. We compare the spectrum 
for the original ensemble {Ihs. plot) to the spectrum for vortex-removed conhgu- 
rations {rhs. plot). Each plot contains the 50 smallest eigenvalues for 10 different 
conhgurations, with every configuration represented by a different symbol. Anti- 
periodic temporal boundary conditions were used for the Dirac operator. 

+ 1 to — 1 when going from a link to its neighbors. Obviously, the spectrum of the 
Dirac operator, which is to a high degree determined by topological properties 
of the gauge field, is sensitive to the discontinuity of the link variables. Thus we 
conclude that analyzing the Dirac spectrum directly for center projected configu- 
rations is a failed, although instructive attempt. We will come back to this issue 
in Section 5. 

3.2 Spectra for vortex- removed configurations 

A way to study the role of center vortices without the drastic measure of project- 
ing the gauge links onto their center elements is the technique of removing the 
vortices outlined in Section 12.21 Although this procedure also introduces short 
range disorder, the changes are by far not as drastic as for center projection (at 
least as long as the vortices are not too fat, see below for more details), and one 
expects the approach via the analysis of Dirac eigenmodes to be viable. 

In the rhs. plot of Fig. |21 we show superimposed the spectra for 10 vortex- 
removed configurations and in the Ihs. plot the spectra for the corresponding 
original configurations. We use the same set of configurations already used in 
Fig. [T] (note that we have changed the vertical scale), and again apply anti- 
periodic temporal boundary conditions. Now the situation in the rhs. plot is 
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quite different from the situation in the rhs. plot of Fig. ^ In particular the 
spectrum does not display the strong shift towards large imaginary parts, and 
essentially stays in the same range as the original spectra. Also the eigenvalues 
are still quite close to the Ginsparg- Wilson circle. However, again we observe a 
gap in the spectrum, indicating that on the vortex-removed configurations the 
chiral condensate vanishes. 

A crucial observation that can be made, is the fact that the eigenvalues fall 
into clusters. It is particularly interesting to count the number of eigenvalues in 
the lower cluster (and its mirror image obtained by reflection at the real axis). 
One finds that this occupation number typically is 8 eigenvalues per configuration 
for this lowest cluster. 

In an attempt to understand this occupation number of the lowest cluster we 
now look at the spectrum for the free case, which can be calculated using Fourier 
transform. Since we are only interested in the lowest eigenvalues, we can ignore 
the effects of the periodicity of the Brillouin zone caused by the lattice and work 
with the continuum form D{p) = ip/^'j^ x I2 of the Dirac operator in momentum 
space. The factor I2 is a 2 x 2 unit matrix coming from the trivial SU(2) field 
configuration. The eigenvalues of D{p) are given by 



Note that each eigenvalue is 4-fold degenerate, where a 2- fold degeneracy follows 
from the block-diagonal structure of the 7-matrices, and another 2-fold degen- 
eracy from the trivial SU(2) color structure. For the momenta we insert the 
discrete momenta allowed on the lattice. These momenta are sensitive to the 
boundary conditions we use. In particular we find (compare Eq. (0) for the 
definition of the boundary condition parameter ip) 



Pi = —ki , i = 1, 2, 3 , P4 = + ip) with = 0, ±1, ±2 ... . (7) 



Here N is the total number of lattice points in one direction and a denotes the 
lattice spacing, i.e. the physical extension L of our lattice is L = aN. 

Let us now analyze what degeneracy of the smallest eigenvalue we find for 
the anti-periodic boundary conditions used in Fig. |21 For this case we have 
ip = 1/2. We find a 4-component of the momentum of p4 = ir/aN for ^4 = 
and p4 = —Tc/aN for ^4 = —1. These two values differ only by a sign, and since 
the eigenvalues ® depend only on the square pi we find an extra degeneracy 
of the lowest eigenvalue (defined by pi = P2 = Ps = 0, p^ = ±TT/aN). Thus for 
anti-periodic boundary conditions we obtain an 8-fold degeneracy as observed in 
the lowest clusters in the rhs. plot of Fig. |21 

It is tempting to interpret the eigenvalues in the rhs. plot of Fig.|21as a slightly 
disturbed free spectrum. If indeed removing the vortices removes all topological 
excitations, then one would expect that the remaining configuration is essentially 
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9 = cp = 1/4 9=1/2 
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Figure 3: Histograms [top row) and cumulated histograms [bottom row) of the 
imaginary parts of Dirac eigenvalues for vortex-removed configurations. We com- 
pare the histograms for periodic [if = 0) b.c. [Ihs. plots), for = 1/4 (center 
plots) and for anti-periodic [if = 1/2) b.c. The vertical dashed lines in the top 
row plots show where we started with summing the numbers for the cumulative 
histograms. In the cumulative histograms [bottom row) we use dotted horizontal 
lines to indicate where the cumulative histogram exhibits a shoulder. 

a trivial configuration plus some fluctuations. In order to test the hypothesis, 
that after removing the vortices one is essentially left with a free spectrum, we 
now use different values for the boundary condition parameter ip. For a value of 
ip = (periodic boundary conditions) we again obtain an 8-fold degeneracy, since 
for Pi = P2 = P3 = Pi = the two signs in © both give A = 0. For the case of 
(/? = 1/4 we obtain only 4-fold degeneracy of the lowest eigenvalue characterized 
by Pi=P2=P3 = 0, Pi = 71/ [2aN). 

To make our study more quantitative we present our results for the spectra 
now as histograms using all 100 conflgurations which we analyzed. In particular 
we use two types of histograms: Histograms where the density of the eigenvalues is 
plotted as a function of ImA, and cumulative histograms where the total number 
of eigenvalues below a certain value of ImA is shown. For the vortex-removed 
conflgurations the results are shown in Fig. El with the regular histograms in the 
top row and the corresponding cumulative histograms in the bottom row. The 
three columns of plots are for the three different boundary conditions we used 
(from left to right: ip = 0,ip = l/4:,(p = 1/2). 

The three regular histograms in the top row are symmetric with respect to 
reflection at the origin due to the corresponding symmetry of the spectrum. The 
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histograms show a clear separation of several maxima, with the position of the 
maxima changing as a function of if. The position of the lowest maxima is 
qualitatively at the same position as for the free case where A = for = 0, 
A = ±i7r/{2aN) for = 1/4 and A = ±i7r/aN for (y9 = 1/2. We remark that the 
50 eigenvalues, available for each of our configurations, are not enough to obtain 
all eigenvalues in the cluster corresponding to the second-smallest free eigenvalue. 
For example in the case of = the second-smallest eigenvalue, characterized 
by one of the being 271 /aN, the others being equal to 0, is already 32-fold 
degenerate. Since this eigenvalue comes as a complex conjugate pair, we find a 
total of 8 + 2 X 32 = 72 eigenvalues in the two smallest eigenvalues (the larger 
one being a complex conjugate pair). Thus with our 50 eigenvalues the second 
cluster is not completely filled and the larger peaks in the top row of histograms 
in Fig. El correspond to only a subset of eigenvalues expected. 

Let us now analyze the occupation numbers in the lowest cluster for the 
different boundary conditions. We start summing our cumulative histograms 
in the minimum below the clusters we are interested in. In the top row plots 
of Fig. El we mark this position by a vertical dashed line. The corresponding 
cumulative histograms are shown right below the regular histograms. For the 
case of V? = we find that the cumulative histogram shows a pronounced shoulder 
at a value of 8 (the cumulative histogram was normalized by the total number 
of configurations). This shows clearly that the lowest cluster of eigenvalues has 
an occupation number of 8, matching the degeneracy of the lowest eigenvalue 
in the free case. For (p = 1/4 the histogram shows shoulders at 4 and 8, again 
matching the degeneracies of the free case. Finally for (y9 = 1/2 we do not see a 
clear shoulder but at least a pronounced dip in the curve. The reason is that for 
this case the lowest and the second eigenvalue are relatively close to each other 
and we do not have a clean separation of the two clusters, as can already be seen 
in the corresponding regular histogram. 

To summarize, the relative position of the peaks in the regular histograms, as 
well as the occupation number in each cluster as obtained from the cumulative 
histograms, leads to the conclusion that the spectrum for vortex-removed config- 
urations strongly resembles a free spectrum, slightly perturbed by fluctuations. 
This finding has two important physical implications: Firstly, we conclude that 
no topological objects are left in the gauge configurations after center vortices 
have been removed. The few zero modes (i.e., real modes on the lattice) that 
are found for ip = (see Fig. ^ are the trivial zero modes present also for the 
free Dirac operator. For the other boundary conditions the complete absence of 
zero modes already shows the absence of topological charge via the index the- 
orem. Secondly, the fact that the spectral density near the origin behaves like 
in the free case shows that the chiral condensate vanishes for vortex-removed 
configurations which strongly supports the findings of OE]. 

For completeness we show the histograms and cumulative histograms also for 
the original configurations (Fig. The regular histograms in the top row do 
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Figure 4: Same as Fig. but now for the original configurations. 

not show any pronounced peaks. The single tall bin at the origin is due to the 
zero modes which are present in the original configurations. The histograms do 
not show any strong dependence on the boundary condition. The summation 
for the cumulative histograms was started at ImA = —0.2 for all boundary con- 
ditions. Furthermore, the cumulative histograms are independent of ip and the 
only feature is the vertical step at ImA = due to the contribution of the zero 
modes. 

3.3 Spectra for random-changed configurations 

In Section 12.31 we have discussed the random-changed configurations which we 
prepared to analyze the effect of an uncorrelated change of the gauge configu- 
ration. For Nfiip = 1000 fiipped links we find that the action goes up twice as 
much as when removing the vortices and it is interesting to see what the effects 
on the spectrum turn out to be. 

In Fig. El we compare the spectra for a random- changed configuration (Ihs. 
plot) to the spectrum for the original configuration (center plot) and the corre- 
sponding vortex- removed configuration (rhs. plot). We used anti-periodic tem- 
poral boundary conditions for all three cases. When comparing the spectra of 
the random- changed and the original configurations, one finds that the former 
is slightly shifted away from the Ginsparg- Wilson circle and slightly compressed 
in the vertical direction. Much more important, however, is the fact that the 
gross features of the two spectra are nearly identical. In particular the number 
of real eigenvalues, i.e. the number of would be zero modes is invariant, and also 
the pattern of the relative spacing of the eigenvalues is very similar. This find- 
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Figure 5: Dirac eigenvalues A in the complex plane. We compare the spectrum for 
the random-changed conhguration [Ihs. plot), the original conhguration {center 
plot) and the vortex-removed configuration {rhs. plot). The random changes, 
as well as the procedure of vortex-removal, both started from the same original 
configuration - the one we also use in the center plot. Anti-periodic temporal 
boundary conditions were implemented. 

ing does not only hold for a single configuration, but also bulk observables such 
as the histograms and cumulative histograms for random- changed configurations 
are very similar to their original counterparts shown in Fig. 0J 

We also analyzed the distribution of the topological charge u in the original 
and the random- changed ensemble. The topological charge u was determined 
from the index theorem using the difference of the numbers of left-handed and 
right handed zero modes. The results are displayed in Fig. El where we show the 
distribution of u for the original configurations in the Ihs. plot and for the random- 
changed configurations in the rhs. plot. Again we observe that the changes due 
the random alterations of the gauge configuration are minimal. This confirms 
an earlier finding [211 where it was demonstrated that the topological features of 
the spectrum of the Wilson Dirac operator are quite stable when adding random 
noise. 

There are certainly many different ways of altering the original configuration 
for comparison of the effects to a removal of center vortices. Such alternative 
prescriptions could e.g. be inspired by a particular picture of QCD vacuum exci- 
tations. We remark that currently we are exploring such alternative modifications 
of the gauge field and their effect on the Dirac eigensystem. The random changes 
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Figure 6: Distribution of the topological charge v. We compare data from the 
original ensemble {Ihs. plot) to the random-changed ensemble {rhs. plot) . 



we use here are merely intended to demonstrate that a non-correlated change, 
which increases the action considerably more than removing the vortices, leaves 
the topological features invariant. 



4 Results for the Dirac eigenmodes 

4.1 Scalar density for zero and near zero modes 

After having explored the spectrum for center projected, vortex-removed and 
random-changed configurations, it is also interesting to study the eigenvectors 
of the Dirac operator for these ensembles. A suitable observable is the scalar 
density p{x) of an eigenvector v. It is defined by summing the color and Dirac 
indices of v at each lattice point x, 

pix) = j:\^ix)j'. (8) 

c,d 

Since the Dirac operator transforms covariantly under a gauge transformation, 
the scalar density p is gauge invariant. 

Since the low-lying Dirac modes do not see large fluctuations of a gauge field, 
the scalar density of the low-lying Dirac modes can be used as a detector of 
(smooth) center vortex flux, which does not rely on topological properties of the 
gauge fields like their topological charge. 

In Fig. [7| we show the scalar density over a 2-d slice through the lattice for 
a zero mode. The Ihs. plot is for the original configuration, the rhs. plot after 
random changes. The slice was chosen such, that it cuts through the maximum 
of p{x). We find a strongly localized lump, which, in an orthodox interpretation. 
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Figure 7: Plots of the scalar density p{x) for the zero-mode of the original con- 
figuration {Ihs. plot) and the corresponding random-changed configuration [rhs. 
plot) . We show the same 2-d slices for both configurations. 

would be the zero mode due to an instanton in the underlying gauge field^ (see 
Refs. j3niEIlE2] for an alternative picture.) The lump is essentially unchanged by 
random changes. Since for the vortex-removed configurations all the zero modes 
are gone, we cannot produce an equivalent plot for this ensemble. 

However, it has long been known that also the "near zero modes" , correspond- 
ing to eigenvalues with small but non-vanishing imaginary parts, show lumpy 
structures^. Since also the vortex-removed configurations do have such near zero 
modes it is possible to compare the scalar density for all four ensembles, original, 
vortex-removed, random-changed and center projected. 

Fig. IHl is such a comparison, with the top left plot showing the original confi- 
guration, the top right plot the random-changed configuration, the bottom left 
plot is for the vortex-removed and the bottom right for the center projected con- 
figuration. Again we show a slice through the maximum of p. The original and 
random-changed configurations show pronounced lumps which are located at the 
same position and have essentially the same shape. When slicing the scalar den- 
sity for the vortex-removed configuration at the same position we find that the 
lump is gone completely. In the plot we even stretched the vertical scale by a 
factor of 10 to make the remaining small wiggles visible at all. Also at other 
positions on the lattice we do not find localized structures in the near zero modes 
of vortex-removed configurations. This holds not only for the particular config- 
uration used in Fig. ISJ but is a generic feature of the whole ensemble. Thus we 
must conclude, that removing the vortices also removes the lumpy structures in 
the eigenmodes. This observation is in agreement with our findings concerning 
the absence of the chiral condensate: The density of near zero eigenvalues which, 
according to the Banks-Casher formula, is necessary to build up the chiral con- 
densate, comes from topological objects which perturb each other only slightly. If 
the chiral condensate is gone for vortex-removed configurations, one expects also 
a dramatic alteration of the near zero modes. This is exactly what we observe. 

^Later we will identify lumps of fractional topological charge with vortex intersection points. 
^According to the instanton picture these lumps originate from instantons and anti- 
instantons perturbing each other. This perturbation is, however, only weak, such that locally 



14 



Figure 8: Plots of the scalar density p{x) for a near-zero-mode of the original 
configuration [top left), as well as the corresponding random- changed {top right), 
vortex-removed {bottom left) and center projected {bottom right) configurations. 
We show the same 2-d slice for all four configurations. The vertical axis for the 
vortex-removed and center projected plot {bottom) is stretched by a factor of 10. 

Surprisingly, in the center projected configuration the locahzed lump is also 
gone. Prom this one might be tempted to conclude that the topological lump seen 
in the original lattice configuration is not related to center vortices at all. How- 
ever, in Section 5 we will provide arguments that, very likely, the toplogical lumps 
are not absent in the center projected lattice configuration but arc just not seen 
by the scalar density of the lowest non-zero mode. The reason is that the center 
projected lattice configurations fiuctuate to wildly to be seen by a low energy 
filter. In Section 5 we will provide arguments supporting this interpretation. 

4.2 Local chirality of near zero modes 

As addressed above, the low-lying near zero modes are expected to locally re- 
semble true zero modes. In particular, they are expected to be locally chiral. In 

order to test for local chirality, one can introduce chiral (right- and left-handed) 
densities p^{x), p^{x), which are obtained by projecting the scalar density p{x) 
to a specific chirality, 

P±{x) = v{x)*^^[l±-f5\dd'v{x)cd' , (9) 

c,d,d' 

where the subscripts denote color and Dirac indices. If the near zero modes are 
locally chiral, then for a given space-time point x only one of the two densities 
p+{x), p-{x) is non-zero. 

the near zero mode still resembles the case of the unperturbed instanton. 
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Figure 9: Distribution of the local chirality variable X. We compare data {from 
left to right) from the original ensemble to data for the random-changed ensemble, 
the vortex-removed and center projected configurations. 



An observable which further analyzes the properties of the near zero modes 
is the local chirality observable X, introduced in [321 Ell studied by several 
groups [331- jHl. The local chirality variable X is obtained by mapping the ratio 

r(x) = ^ , (10) 

which takes values in the interval [0, oo), to the interval [—1, 1], 

X{x) = ^ tan^i (v^K^) - 1 • (11) 

If the near zero modes are locally chiral, one expects that the distribution of X 
shows a double peak structure with maxima near ±1. If on the other hand the 
near zero modes are not locally chiral, one expects a single peak near the origin. 
We determine this distribution using histograms for the values of X{x), where 
we include all those lattice points x where we find the largest 12.5% of the scalar 
density. Thus we analyze the local density only for the highest peaks in p{x). 

The results for the local chirality are shown in Fig. [HI For the original and 
the random-changed configurations, we find a clear double-peak structure indi- 
cating that the near zero modes are locally chiral for these two cases. For the 
vortex-removed configurations the double-peak structure is gone, indicating that 
the wiggles in p{x) that remain after removing the vortices are not locally chiral 
structures. This confirms our previous interpretation of these wiggles as struc- 
tureless fluctuations. As before, we find that the center projected configurations 
are too singular for an analysis with the low-lying Dirac eigenmodes. 
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5 Results for smooth center vortices 



During the course of our analysis we have found that center projected config- 
urations do not give sensible results when analyzed with the low-lying Dirac 
eigenmodes. We have argued that center projection yields configurations that 
are too singular for the eigenmode analysis. This is a consequence of the center 
projection, which converts the originally fat center vortices (present in the full 
Yang-Mills ensemble) into "ideal" center vortices, whose transversal extension is 
one lattice spacing only. 

However, in this Section we demonstrate that sufficently smooth center vor- 
tices do indeed produce zero modes and thus can also give rise to chiral symmetry 
breaking. For this study we use smooth vortex configurations consisting of two 
pairs of parallel smeared out planar vortex sheets, intersecting perpendicularly 
in four (smeared out) intersection points, each carrying the topological charge 
±^ (for more details see Ref. [lEl)- The vortex sheets are closed by the periodic 
boundary condition on the 4-torus. The orientation of the vortex sheets, i.e. the 
direction of the flux is such that the contributions from the intersection points 
add up to a total topological charge of u = 2 for configuration #0 and i/ = for 
configurations #1, #2, #3, #4, respectively. 

5.1 Spectra of smooth center vortex configurations vs. 
center projected vortices 

Let us begin with analyzing the Dirac spectra for our smooth center vortex con- 
figurations. In the Ihs. plot of Fig. we show the Dirac spectra for our 5 
configurations and compare them to the spectra obtained for the corresponding 
center projected configuratons. 

The Dirac spectra of the smeared out center vortex configuration have no 
gap around zero virtuality. In particular, for configuration 7^0 with non-zero 
topological charge u = 2 there is a zero mode, which is two-fold degenerate. After 
center projection the gap in the Dirac spectra emerges and the zero mode is gone. 
This is not surprising since in the process of center projection the center vortices 
lose their orientation and thus their total topological charge, although they still 
carry the individual local spots of fractional topological charge. But these local 
contributions will usually add up to zero total charge. This is because orientable 
closed surfaces in D = 4 have zero total self-intersection number, which up to a 
factor of I represents the topological charge of center vortex surfaces [THl IT^. 

To illustrate the effect of center projection and center vortex removal in more 
detail we compare in Fig. ^2 the Dirac spectrum of the fat center vortex configu- 
ration #0 (see Table ^ with topological charge u = 2 with those of its (b) center 
projected and (c) center vortex-removed counter parts. For sake of comparison, in 
Fig. ITTT d) we also show the Dirac spectrum of the fat center vortex configuration 
#4, which has the same vortex surfaces as the configuration #0 (a), however. 
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Figure 10: Dirac spectrum for smooth center vortices (Ihs.) and their center 
projected counter parts {rhs.). Note that the center projection is the same for 
configurations #0 and #4. 



with different orientations of their ffuxes such that z/ = (cf. Table Q). Both 
fat center vortex configurations #0 and #4 have the same center projected and 
center vortex-removed counter parts and thus the same Dirac spectrum shown in 
Fig. lTlT b) and (c). As is seen, center projection (b) changes the Dirac spectrum as 
drastically as center vortex removal (c). This is somewhat counter intuitive since 
the configurations under consideration are plain vortex configurations but, in fact, 
can be easily understood: The vortex configuration considered in Fig. ITTT a). (d) 
are extremely fat. In the transversal directions they stretch out over the whole 
lattice. Accordingly, there are only small gradients involved and our lattice Dirac 
operator works well: all eigenvalues are on the Ginsparg- Wilson circle, see fig 
ITTT a). (d). This property is lost in the center projection, which replaces the 
originally fat vortices by the ideal (very thin) ones, whose transversal extension 
is one lattice spacing. Hence the center projected configuration contains large 
gradients which cannot be properly captured by our Dirac operator. Accordingly 
its eigenvalues are off the Ginsparg- Wilson circle. The method of center vortex 
"removal" subtracts the ideal center vortices arising in the center projection from 
the fat vortices. This will give rise to even larger gradients and accordingly the 
Dirac spectrum Fig. ITlT c) will show more "noise." The method of center vortex 
removal obviously works the better the thinner the vortices, and becomes perfect 
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conf. 


vortex 1 


vortex 2 


vortex 3 


vortex 4 


u 




(^0, h, f) 


(^0, h, f) 


(^1,^2,/) 


(^1,^2,/) 




#0 


(10,14,0.5) 


(4, 14,0.5) 


(4,4,0.5) 


(12,14,0.5) 


2 


#4 


(10,14,-0.5) 


(4,14,0.5) 


(4,4,0.5) 


(12,14,-0.5) 






Table 1 : Smooth center vortex configurations consisting of 4 planar center vortex 
stieets wfiich are all parallel to coordinate planes ki — kj. The table shows the 
flux f and the coordinates (fco, k^) of the two vortices parallel to the ki — k2 plane 
and the coordinates {ki, ^2) for the two vortices parallel to the k^ — k^ plane for 
the conhgurations #0 and #4. Conhguration 7^4 arises from conhguration ^0 by 
changing the orientation of the flux. Furthermore, the total topological charge v 
of the configuration is given. 



for ideal center vortices (having a transversal extension of one lattice spacing)^. 
The thinner the vortex the less flux is left after vortex removal. To illustrate this 
we show in Fig. the Dirac spectra for the vortex configurations which result 
from those considered in Fig. II II when their transversal vortex extension is shrunk 
to two lattice spacings. The starting center vortices now have larger gradients 
and the Dirac eigenvalues are now shifted off the Ginsparg- Wilson circle. At the 
same time the Dirac eigenvalues of the vortex-removed configuration Fig. (c) 
are now basically on the Ginsparg- Wilson circle implying that this configuration 
contains only little gradients contrary to what we have observed for the very fat 
vortices (Fig. ^2(c)). For sake of completeness let us also mention that the cen- 
ter projected configurations of the very fat and thin center vortices are the same 
and so are the corresponding Dirac spectra shown in Fig. ^2 (b) and Fig. (b), 
respectively. 

5.2 Eigenvectors 

After having demonstrated that for smooth center vortices the spectrum indeed 
matches the expectation from the index theorem, let us now analyze the corre- 
sponding eigenvectors. 

Fig. ^1 shows the scalar density of the lowest lying non-zero mode of the fat 
vortex configurations #0 {u = 2) and #4 {u = 0) and their center projected 
image. Both configurations have up to orientation the same same vortex flux 
and thus the same center projected image. The scalar density is shown for a 
plane which does not coincide with one of the vortex planes. As is seen the 
scalar density of the lowest non-zero mode detects only one of the two parallel 
ki — k2 flux sheets at fco = 4 and ko = 1, and furthermore different sheets for the 
u = 2 and u = configurations. The scalar density of the lowest Dirac mode 

•^The method of vortex removal converts an ideal center vortex into the trivial configuration 
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Figure 11: Dirac spectra of {a) the fat center vortex configuration ^0 {u = 2), its 
[b) center projected and (c) vortex-removed counter part, and [d) the fat center 
vortex conhguration with u = obtained from (a) by changing the orientation of 
the flux. 
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Figure 12: Dirac spectra for the same fat vortex conhguration as in FigMBexcept 
that the transversal extension of the vortex flux was shrunk to 2 lattice spacings. 



of the center projected configuration sees both center vortex fluxes and in more 
pronounced form than the corresponding Dirac mode in the original fat center 
vortex background. Furthermore, the scalar density is concentrated at the vortex 
intersection points. These results are in accordance with the findings of Ref. [Tl)] 
where the concentration of the scalar density of the zero modes near the vortex 
flux and, in particular, at the intersection points was found. Surprisingly, the 
scalar density of the center vortex-removed configuration is also localized at the 
vortex flux and is basically the same as for the center projected configuration. One 
would have expected that the localized structures arising from the center vortices 
and their intersections are eliminated by the vortex removal. However, one should 
keep in mind that vortex removal by the method of Ref. [3], implies subtraction 
of an ideal center vortex from a fat vortex in the original configuration. If the 
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original vortex is very fat, vortex removal near the vortex core is almost the same 
as center projection. Only for thin smooth center vortices the vortex removal 
procedure properly removes the vortex flux. 

Fig. El shows the scalar density for the same configurations as in Fig. IT!?| 
however, for a different 2-d slice which coincides with a vortex plane. In addition, 
the transversal size of the vortex flux was shrunk to two lattice spacings. Again, 
the scalar density is localized at the intersection points and is basically the same 
for original and the center projected configurations while it is supressed for the 
vortex-removed configuration, in particular at the intersection points (Note the 
different scales!). 

Thus for very smooth fat center vortex configurations the center projected 
vortex surfaces are also smooth and their fluxes and topological spots are very 
well seen by the scalar density of the near zero modes. 

Contrary to the scalar density p{x) (jSj), the chiral densities p±{x) © do 
feel the orientation of the flux. This can be seen in Fig. E[ where we show 
the chiral densities of the lowest non-zero Dirac modes of the configurations 
#0 and #4, which differ only in the orientation of the vortex flux. The chiral 
densities are shown in a 2-d slice in which these two configurations have two vortex 
intersection points. For configuration 7^0 the two vortex intersection points have 
the same orientation and accordingly give both rise to peaks of the same size in 
the chiral densities p±(x). However, the density p+(x) is suppressed by a factor 
of 20 compared to p_ (x) implying that the considered non-zero Dirac mode is at 
both vortex intersection points approximately left handed (negative chirality) in 
agreement with the vortex orientation. For configuration #4 the two intersection 
points have opposite orientation and we observe a peak in either p+(x) or p_(x) 
at the two intersection points. The orientation of the flux is lost in the center 
projection which hence converts both center vortex configurations #0 and #4 
into the same ideal center vortex. Therefore the center projected configurations 
do not distinguish between p+{x) and p-{x) and consequently do not reproduce 
the chiral densities properly. This is indeed seen in Fig. where we compare 
the chiral densities of the lowest lying non-zero Dirac modes for the fat center 
vortex configurations 7^0 and ^4, the corresponding center projected and the 
center vortex-removed configurations. 

6 Summary and Conclusions 

Using the chirally improved Dirac operator constructed in [2011^; we have stud- 
ied the influence of center vortices on various properties of the Dirac spectrum. 
We have shown that removal of center vortices eliminates the zero modes and near 
zero modes of the Dirac operator. Via the Banks-Casher relation this implies the 
restoration of chiral symmetry. We have demonstrated that the spectrum for 
vortex-removed configurations strongly resembles the free Dirac spectrum. For 
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Figure 13: Plots of the scalar density of the lowest non-zero eigenmodes of con- 
Ggurations #0 {upper left) and #4 {upper right) and of their center projected 
{lower left) and center vortex-removed counter part {lower right) in the 2-d slice 
defined by {ki,ks) = (12,14). In this plane the configuration has intersection 
points at {ko, k2) = (4, 14) and (10, 14). 




Figure 14: Plots of the scalar density as in Fig. [131 but in the 2-d slice de- 
hned by the second vortex, i.e. {k^^k^) = (4,14), see Tabled In this plane 
the configuration has vortex intersection points at {ki, /C2) = (4,4) and (12, 14). 
Furthermore, we show the scalar density of the first non-zero mode for the corre- 
sponding vortex-removed conhguration {lower right) and of the center projected 
conhguration {lower left) . 
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Figure 15: (a) Plots of the chiral densities p+{x) {upper row) and p-{x) {lower 
row) of the lowest non-zero Dirac mode for the configuration #0 and its center 
projected and vortex-removed counter parts through the 2-d slice defined by the 

plane {ko,ks) = (4,14), in which this configuration has two vortex intersection 
points. Note that in the upper left plot the density was stretched by a 

factor of 20. (b) The same as in (a) for conhguration ^4. 
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the Dirac eigenmodes we found that removing the center vortices destroys local 
chirality of the low lying modes. Comparison with configurations, subject to 
incoherent random changes, shows that removing the center vortices specifically 
targets the topological properties, while random changes leave them intact. 

Somewhat surprisingly we find that also center projection eliminates the topo- 
logical charge and gives rise to a gap near zero virtuality. However, this result 
does not mean that center vortices are not responsible for the topological proper- 
ties of the gauge fields and for the spontaneous breaking of chiral symmetry but 
can be traced back to the fact that center projection removes the orientation of 
the vortices. Center projection keeps track, however, of the module of the flux 
of the vortices. This was demonstrated by measuring the scalar density of the 
lowest-lying non-zero modes, which is insensitive to the direction of the flux. The 
center projected configurations show the same localization of the scalar density 
near the center vortices (and their intersection points) as the full configurations 
when the latter are rather smooth. While the center projected configurations 
reproduce the localization of the scalar density near the vortices, they do not 
reproduce the localization of the chiral densities, as the latter depend on the ori- 
entation of the flux. The orientation of the flux of the center vortices is irrelevant 
for their confining properties but crucial for their topological and chiral prop- 
erties. Therefore we arive at the conclusion that the familiar method of center 
projection is not suited to study the chiral and topological properties of center 
vortices for lattice ensembles. To capture the topological properties one has, at 
least, to embed the center configurations into the Cartan subgroup to preserve 
the orientation of the center surfaces. 
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